ee PUBLISHED BY THE UNIVERSITY OF ILLINOIS 
RPE vice as URBANA, ILLINOIS 


74 : 
an Team 
ee Pek \ fae ; 
Rees y Le i 4 L 
Ltd 1? te . 


Bs {HE Engineering Experiment. Station was cétablished by act a 
je | ia of the Board of Trustees of the University of Illinois on De- « 


cember 8, 1903, It is the purpose of the Station to conduct 


ae ? investigations and make studies of importance to the engineering, M 
manufacturing, railway, mining, and other industrial interests of the ~ 


State. | 
_ The management of the Engineering. Experiment Siation is vested 


in an Executive Staff composed of the Director and his Assistant, the e 


Heads of the several Departments in the College of Engineering, and 


the Professor of Chemical Engineering. This Staff is responsible for ; 
the establishment of general policies governing the work of the Station, ~ 


including the approval of material for publication. All members of 
‘the teaching staff of the College are encouraged to engage in scientific 


research, either directly or in codperation with the Research Corps, 


composed of full-time research assistants, research graduate assistants, 
and special investigators. 

To render the results’ of.its scientific investigations available ia 
the public, the Engineering Experiment Station publishes and dis- 
tributes a series of bulletins. Occasionally it publishes circulars of 
timely interest, presenting information of importance, compiled from 


various sources which may not readily be accessible to the clientele, 


of the Station, and reprints of articles appearing in the technical press 
written:by members of. the staff and others. 


The volume and number at the top of the front cover page are | 


merely arbitrary numbers and refer to the general publications of the 


University, At the top of the inner title page is given the number of. - 


the Engineering Experiment Station bulletin, circular, or reprint which . 


should be used in referring to these publications. 
For copies of publications or for other information address 
THE ENGINEERING EXPERIMENT STATION, 
UNIVERSITY OF ILLINOIS, 


Urzana, InLinots 


a eae 
hee 


- 


Burren No. 309 


- 


IN THE ARBITRARY CONSTANTS O 
SIMULTANEOUS EQUATIONS 


BY 


GEORGE H. DELL 


ASSOCIATE IN CrvIL ENGINEERING 


a gs 


PUBLISHED BY THE UNIVERSITY OF ILLINOIS 
ENGINEERING EXPERIMENT STATION 


Price: Srxry Cents 


a 


3000—12-38—15611 : ' 4 guversiry ; 
it PRESS 


: 2. Si ieavwiadenenet 


a 3. Previous Contributions . 
oo of the Investigation 
na i 7 
Be Basic : ConcEPTS 
ms, Definitions 
. 6. General Relations 


II i PRIMARY ERRORS AND PROBABLE ERRORS . 
a 
8. Solution by Determinants . 
Example’ No. 1 
9. Solution by Elimination 
Example No.2. 
10. Solution by Means of Unit “tee iss = at 
Example No. 3 oe ei A 
11. Finding the Maximum eee. eax 
Example No. 4 ra 
12. Primary Errors in the etnowis of otiincars 
Equations 
Example No.5. 
= 13. Effect of Errors in Right ad Tokuias 
Example No. 6 ea 3 : 
14. Partial Differentials and Par Hal Detvatives ; 
15. Probable Errors . 
Example No. 7 
Example No. 8 


IV. True Errors . : 
16. Methods ible fer Cleatine re gene abe . 


17. The True Maximum Error . 
Example No. 9 
: Example No. 10 
| 


ye Primary Errors in the Unknowns of Linear Paustions 


c~. . | 
_ VI. Conctustons ae 
---:19, Comparison of Results’ 
20. Conclusions . . . 


APPENDIX Ay (i) *r-) 4g ot OR 2 
1. Errors of Higher Orders in the Unknowns of __ 
Linear Equations . . . . . . 


THE EFFECTS OF ERRORS OR VARIATIONS IN 
THE ARBITRARY CONSTANTS OF 
SIMULTANEOUS EQUATIONS 


I. IntTRopDUCTION 


_ 1, Introductory—Many problems in science and engineering in- 
volve the use of quantities determined by means of simultaneous 
equations whose coefficients are subject to errors or variations. 


Examples in the field of structural engineering occur in the analysis 


of continuous frames, where the coefficients of the governing equa- 
tions depend upon various ideal dimensions and elastic properties 


which are not exactly realized in the actual structure. A similar 


situation exists in connection with the determination or the interpre- 
tation of experimental data, if equations which are subject to observa- 


tional errors or physical variations are involved. It is appropriate, 


therefore, that an inquiry be made into the effects of such errors. 

The basis of the problem consists in the presence of errors of 
unknown signs and magnitudes in the arbitrary constants entering 
into the equations. These arbitrary constants are generally of a 
geometrical or physical nature, and include parametric, dimensional, 
and empirical constants. 

In studying the effects of such errors upon the computed values 
of the unknowns, or upon given functions of the unknowns, the aim 


of the investigation may be either to determine the probable errors, 


or to calculate the maximum effects which may be produced by errors 
of assigned magnitudes in the arbitrary constants. In the former 


ease, the probable errors of the arbitrary constants may be estimated 


or obtained by experiment; in the latter case, the magnitudes to be 
assigned to the errors in the arbitrary constants may be determined 


with reference to the corresponding probable errors, or may be 
_ assumed arbitrarily. 
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the material. 


3. Previous Contributions —Professor Blumenthal, at Aachen,* in 


investigating the precision of the solution of linear equations, derived 
an expression by means of which he was able to determine an upper 
bound to the error in the unknown which was most seriously affected 
by errors of restricted magnitudes in the coefficients. Professor 
A. Hertwig, also at Aachen,{ derived an expression for the primary 
error in a given unknown of a system of linear equations in which 
the errors in the coefficients were taken at a fixed percentage, and 
showed how to determine the signs in order to obtain the maximum 
primary error. He also stated the fundamental formula defining the 
mean error of a given unknown (the mean error being a quantity 
related to the probable error) and derived an approximate method 
for evaluating the same. In addition, he investigated the ‘‘error- 


sensitiveness”’ of a given system of equations, applied to the analysis 


of statically indeterminate structures, with a view to determining 
the relative advantages of choosing different quantities as the redun- 
dants in the solution. For this purpose he used as a criterion the 
error in the determinant of the coefficients of the governing equations. 


4. Extent of the Investigation—The present treatment is believed 
to include several new features, among them being the derivation of 
general expressions for the primary error by the method of total 
differentials, the calculation of errors of higher orders, the distinction 
between independent and interdependent coefficients, and a discus- 
sion of the true maximum error. The investigation is also extended 
to non-linear equations. 

The study deals specifically with systems containing n linearly 
independent equations in n unknowns, as follows: 

(a) Linear equations with linear arbitrary constants. 

(b) Non-linear equations with non-linear arbitrary constants. 
No difficulty should be experienced in extending the treatment to 
linear equations with non-linear arbitrary constants or to other 
combinations. 

The investigation is particularly concerned with the calculation 
of maximum errors, but since, in many applications, the probable 


f 


*Blumenthal, Ueber die G igkeit der W i i i 
Pie, Base tee enauigkeit der Wurzeln linearer Gleichungen, Zeitsch. f. Math. u. 


: - Hertwig, Die Fehlerwirkungen beim Aufloesen lin Glei i 
statisch unbestimmte Gebilde, Der Hisenbau, Vol. 8, 1917, pp. 110-132, Mang datret meget OE: 
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er ‘ors may be of primary interest, and since a considerable portion 
| of the procedure is common to both quests, the steps involved in the 
determination of the probable errors are also briefly outlined and 


: illustrated. 


II. Bastc Concepts 


5. Definitions—For purposes of definition, consider the following 
system of linear equations: 


fi = Quti +. 2... + dimtm +.... + Ointn — ws = 0 (1) 
fm = Gmit1 +... . + Ommim +. ... + Anatn — Lm = 0 (1m) 


Sn = Gnitr +... + Onmim +... + Anntn — pn = 0 (1n) 


where the coefficients ay .... Gan contain the arbitrary constants 
k, .... kp, which are subject to the errors or variations represented 
by the increments Ak, .... Ak,, respectively; and y1.... un are con- 
stants usually appearing on the right-hand side of the equations, and 
which hence will be referred to as the “right-hand terms.” For the 
present, these quantities will be regarded as absolute constants. 

_ Independent Coefficients. If none of the arbitrary constants 
ky .... kp appears in more than one of the aggregate number of the 
coefficients in the given system, the coefficients are defined as being 
independent. 

_ Interdependent Coefficients. If any or all of the arbitrary constants 
k, .... k» appears in two or more of the aggregate number of the 
coefficients in the given system, the coefficients in question are defined 
as being interdependent. 

Initial Values. The initial value of a given quantity is that value 
determined by using the nominal or mean values of the arbitrary 
constants. The initial values of the arbitrary constants, of the 
coefficients, and of the unknowns are denoted by ki.... kp, ai.... 
Ban, and %;.... %,, respectively. 

Modified Values. The modified value of a given quantity is that 
value determined by assigning definite signs and magnitudes to the 
errors in the arbitrary constants. The modified values of the arbitrary 
constants, of the coefficients, and of the unknowns are denoted by 
Me Gi... - Onn, aNd 21... . Zp, respectively. 

: Original Equations. The original equations are those in which the 


i s 


a, 
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coefficients have their initial values, ati... . @nn; their solution deter-— 
mines the initial values of the unknowns, z?.... x2. a . 
Revised Equations. The revised equations are those in which the 
coefficients have their modified values, ai .... Grn; their solution 
determines the modified values of the unknowns, zi.... 2). a 
True Errors. The true error in a given quantity is defined as the 
difference between the modified and the initial values of that quantity. 
The true error in a given coefficient, a%;, is 


° , o, 
Aai; = ij — aii; 


similarly, the true error in a given unknown, 7%, is 
Aas, = 2), "— 22, 
The true error may generally be represented by a Taylor’s series. 
Primary Errors. The primary error in a given quantity is the 
total first-order differential of that quantity. It is equivalent to the 
first term of the Taylor’s series which expresses the true error. The 
symbol 6x, is used to denote the primary error in z°,. 
The primary error is expressible as a linear function of the errors 
Ak, .... Ak» with the aid of the initial values, k2....k°, a%.... Gnaa 
SNC ice eels 


: 


Maximum Primary Error. The maximum primary error in eo 
given quantity is the numerical value of the primary error resulting | 
from assigning to the errors in the arbitrary constants such signs as | 
will make the separate terms in the expression for the primary error | 
cumulative in their effects. . 

(Note: The term ‘‘maximum” is used here and in subsequent | 
definitions to denote “greatest,” though not necessarily in the sense 
of an extremum, that is, where the derivatives vanish. ) 4 


Errors of Higher Orders. Following the primary error, the succes- | 
sive terms in the Taylor’s series for the true error in a given quantity 
represent the errors of higher orders. The symbols tm, O°, CLC 
are used to denote the errors of higher orders in ope | 

Apparent Maximum Error. The apparent maximum error in @ 
given quantity is the true error obtained by using with the errors in | 
the arbitrary constants a particular combination of signs which is | 
presumed to yield the maximum error. : 

True Maximum Error. The true maximum error in a given quan- 
tity is the true error obtained by using that particular combination 
of signs with the errors (of limited magnitudes) in the arbitrary con- 
stants which will yield the greatest possible error in that quantity. 


d Brror-coeficients. The coefficients of the increments Bhy ehh; 
he expression for the primary error of a given quantity, ated in 
nerical or in literal terms, are called “error-coefficients.” Since 


Benaly error is identical with the total first-order differential, 


EProbable Error. The probable error is a Pikaties pertaining ‘ 
the probability curve or error function. It is defined as that error 
which is Just as likely as not to be exceeded.* 


4 6. General Relations.—If from a function z(r, s, t,....) of several 
ir dependent variables r, s, t,...., itis desired to develop the function 


x + Ax® = 2(r° + Ar, 8° + As, tf? + At,....), 


where 2° = x(r°, s°, t°, ....), the indicated expansion is given by the 
ollowing Taylor’s series: 


Ox Ox Ox 
a? + Ax? = 4 (Sar¢ as eee) 
or Os ot 


1 /dz @) 
rel = as ne Me 
Dee ge “s + 
1 /dz (n) 
ae 14Y OF (n) 
or, ia art as+ art. y +R, (2) 
;- n=1 n! or ° 


is used to express sym- 


°o 


Ox Ox Ox (n) 
where (ar + — As sewage +... .) 


or Os 


‘bolically the n-th order total differential of z, evaluated by means of 
‘the BTIIALVAIUCS,, 105189) 1, 6 5 aes t 

It is assumed that z is continuous and that all of the derivatives 
ised exist and are continuous. If zx has a finite number, n, of deriva- 
tives, Rn, = 0 (Taylor’s Formula). If there are an infinite number 
of derivatives it is assumed that R, approaches zero, and that the 
‘series is absolutely convergent throughout the region considered 
(Taylor’ s Series). 


*See Art. 132, Higher Mathematics for Engineers and Physicists, by I. S. and E. 8. Sokolnikoff, 
fArt. 51, Goursat-Hedrick, , Math. thal Vol. I; oe 2, Chap. VII, Osgood, Advanced Calculus, 
tArt. 15, Goursat-Hedrick, Math. Analysis, Vol. I 
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a 
Equation (2) may also be written in the following forms: E 
ax d"x my 3 
Ax® = dz} + +eeeet ee dnp ay (3a) 

° 2! ° n! ° ‘ 
and Axe = 62+ 8a+....+6%4+.... (3b) 


where, by definition, the primary error is 6x = dz 


, the secondary 
Cx ‘ 


oO 


error is &2 = ,» ete. 


oO 


Numerical examples show that Equations (3a) and (3b) may be 
extended to composite functions, of the type x(r, s,t,....5u,0,....), 
where r,s, ¢,... . are functions of the independent variables u, v,..... 


III. Primary Errors AND PROBABLE ERRORS 


7. Primary Errors in the Unknowns of Linear Equations.—Let 
Equations (1) to (1n), inclusive, in Section 5, represent a system of n 
linearly independent linear equations in n unknowns, where the 
coefficients ai... .. nn are integral linear functions of one or more of 
the p arbitrary constants k, .... k, occurring in the given system, 
which are subject to the errors or variations Ak, . .. . Ak», respectively. 

A typical coefficient may be represented by the function aij(ka, 

., Ky), whence 

003; 04;; 
Aa;; = da;; = Ak #....4+ 


°o 1 Dp 


Ak, (4) 
It is assumed that the magnitudes of the errors in the arbitrary 
constants are so limited that the equations will remain linearly inde- 
pendent (that is, the determinant of the coefficients does not vanish), 
and that the Taylor’s series representing the true error in any of 
the"unknowns will be absolutely convergent. , 
te} From Equations (1) to (1n), inclusive, the total differentials dfi, 
erdfad. Ach Ufo areas. thilowse 


df, = anda, + Aas 8 + AimALm ae ei na} + Andy 
“b tidon +... 22 ada, +. dae 0 (5) 


Om = Omidt, +... + mmm. ... + nna 
+ 21dOmi +... : + tmddmm +... + taddnn = 0 (5m) 


ERRORS IN SIMULTANEOUS EQUATIONS. 11 


df, = Gudz+....+ fase peek + Onndtn 
oo Xidani sw <tr LmdAnm +. AG Spies =0 (5n) 


‘From these equations the corresponding relations between the errors 


in the coefficients and the primary errors in the unknowns may be 


obtained by substitution of the initial te Gini ieee Genaand 


Hi.... 23, as follows: 


aj,éa, +....+ a?,,62,, - .... +252, + < x°Aa,; = 0 (6) 
jaa 

O5100, + .... + ay,ot,, +.... 7 42,02, + >) a%dc,,; = 0 (6m) 
j=1 


Str)... 4 o,f, +.... + 0°,02, 42> 2%Ac,,=0 (6n) 
j= 


In solving for the primary errors 621 .... 6%n, it is desirable, in 
the first stage, to employ an abbreviated symbol for the terms indi- 
cated by the summation sign in the foregoing equations. For this 
purpose, let 

= dX x°Aa;; (7) 
ry ee 
The following equations may then be written in place of Equations 
(6) to (6n): 


ST er ee See a ee (8) 


O91601 +... + Oembim +... + Omndtn = —Cm (8m) 
07162; +... . t+ Onmdtm t+... . + Anndtn = —Cn (8n) 


It will be noted that the left-hand sides of these equations differ 


from the original equations only in containing 62, .... 5x, in place 


of the unknowns 2 ....%n; and that the right-hand terms represent 
the differential changes produced by the errors or variations in the 
arbitrary constants of the given equations, multiplied by —1. 

A solution of Equations (8) to (8n) will give the desired expressions 
for the primary errors 521 .... 6%n, in terms of @.... @n. Various 
methods of obtaining these expressions will be discussed and illus- 
trated in the following sections, before proceeding with the determina- 


tion of maximum errors. 


’ © 
t ? , 
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inclusive, Section 5, the determinant of the coefficients, D, is as 
follows: RG 


~. 


Cyr tata a taiaheee Oe 
POSE at Gee (9a) 


oO fo) fe) 
Gist ssh Coane ce ee 


_ The initial value of a given unknown, zm, may be obtained from — 
the relation . 
. In 


ae, (10a) 


Oy ax lige se eles 


Noting now that the determinant of the coefficients of the é2’s in _ 


a Equations (8) to (8n), inclusive, is identical with D, as given by : 
we Equation (9a), the following solution is obtained for 62m: p 
} Oia ss Clips nae, } 
5 : | (11a) 
Im = ini wile Og nein «hee 
Sine 
OE 2 Rs 


mn 


A comparison of the expression for 2°, Equation (10a), and the 
corresponding expression for 62m, Equation (11a), shows that the 
expression for 62» is obtainable by substituting the e terms in place : 
of the uw terms in the former expression and dividing by —1. This — 
suggests the development of the determinants indicated by Equa- 3: 
tion (10a) with respect to the elements of the m-th column. In order — 
to avoid errors in sign in transforming the expression for zs, into that 
for 52m, it is advisable to preserve the signs of the » terms as they 
appear on the right-hand side of the equations. This procedure is 
facilitated by the use of signed minors, or “‘co-factors.”* For example, 
in Equations (1) to (1n), inclusive, the coefficient in the 7-th row and 


*Art. 17, Higher Mathematics for Engineers and Physicists, by I. S. and E. S. Sokolnikoff, 
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m-th ¢ olumn o ihe determinant D is a%m, God its re minor 


“5 ) i+ “Aim = = Can 
Ji he expressions for 2 e. and 62m May then be written as follows: 


_ 


D = oy BimC im (9b) 
t=1 ’ 
. 4 = Dp pe (10b) 
and bXm eo Wee yo re Cen i ' 3 J 
ay —> e; (11b) 


Example No. 1. It is desired to investigate the effects of errors 
in the coefficients of the following system of equations, of the type 
Kkax + ky + kise = wi: 


22° + 3y° + 32° = 12.85 (12a) 
x° — 5y° + 22° = —13.10 (12b) 
e+Qaye—- 2= 8.05 (12¢) 


Using the method of determinants, obtain the expression for the 
_ primary errors 62, 6y, and 6z, in terms of é, @, and és. 

Solution: In accordance with Section 7, the following equations 
may be written: 


25x + 36y + 362 = —e (13a) 
az 6x — 5dy + 262 = —e (13b) 
6a + 2Q5y — 62 = —@ (13¢) 
where é: = t°Aky + y°Aki2 + 2°Akis 
€. = £°Aka: + y°Akog + 2°Akos 
and €3 = r°Aks, + y°Aks, + 2°Aks33 


The determinant of the coefficients of Equations (12) and (13) is 
7 Es: 
Des 15 2 
1 2. —1 


4 1 °en Je 


- 


BOS ing, Raf as ae y! fas pyaats ; 
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The development of this determinant with respect to the eleme nts 


of the first, second, and third columns, respectively, by means | of 
co-factors, is as follows: — ' 


Dr= 20) #1C 9)-410 21) = 32 
Dz = 3(3) — 5(—5) + 2(— 1) = 32 
D; = 3(7) + 2(—1) — 1(—18) = 32 7. 


the co-factors being shown in parentheses. “A 
In accordance with Equations (9b) and (10b), the unknowns of | 
Equations (12) are 4 


~ 


1 
: = Bq [12-85(1) — 13.10( 9) + 8.05( 21)] = 2.00 


aa [12.85(3) — 13.10(—5) + 8.05(— 1)] = 3.00 


1 
and 2° = a [12.85(7) — 13.10(—1) + 8.05(—13)] = —0.050 


The desired expressions for the primary errors in the unknowns, 
as given by Equation (11b), are then as follows: 


7 ; 
‘a . ove 39 ( ey + 9e, + 21es) 


1 
jy = ae (3e, — 5e2 — es) 


and (a= 


39 (Te, —— ae 13¢3) 
which is likewise the solution of Equations (13). : 

9. Solution by Elimination.—In case it is practicable to solve the ¥ 
original equations by the method of elimination, a parallel solution, a 
in which the right-hand terms of the original equations are replaced 
DY p=" 6124 en eral give the expression for the primary errors, = 
621... . 5 (see Equations (8)). & 


Example No. 2. Derive the expression for the primary errors of - : 
the unknowns of Example No. 1, in terms of 1, e, and es, by solving — 
Equations (13) by the method of elimination: : 


> 
a 


a Sb aes ea ah Ce RO Re ao et SORTS ee alr 
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(15 


26x + 35y + 352 = —e, (18a) 
6x — 5dy + 26z = —e (13b) 
dx + 2y — 62 = —es (13¢) 


Solution: Elimination of the first unknown from Equations (18a) 
a nd (13b), and from Equations (13b) and (13c), results in the 
equations 


13éy — 62 = —e, + 2e (14a) 7 
p —Toy + 362 = —e+ 5 (14b) 


Elimination of 6z from these equations gives 
i 
by = 32 (—3e a 5é2 + €3) 
Substitution of this value in Equation (14a) gives 
1 
6z = 1386y +e: — 2 = ee (—7e + e + 1363) 


Finally, by substituting the foregoing values for é6y and 6z in Equa- 
tion (13b), one finds 


1 
oC = =x a 5by — 262 = 32. (-—e a Jen — 21es) 


10. Solution by Means of Unit System (um = —1).—The coeffi- 
cients of e,.... é, in the expression for the primary error 6z,, in a 
given unknown z, may all be obtained by solving a new system of 
- equations, in which the right-hand term of the m-th equation is 
_ placed equal to —1, and the corresponding terms of the remaining 

equations are zero. For convenience, such a system will be referred 
to as a “unit system.” 
The left-hand sides of the equations in the unit system are formed 
by replacing the elements of each row in the determinant of the 
- coefficients of the original equations by the elements of the corre- 
sponding column. 

By way of demonstrating this method, consider the following 

_ system of equations in three unknowns: 


A121 + Aree + A303 = m1 
Qgi%1 + AgeX2 + A23%3 = me 


3101 + AgeX2 + A333 = Ms 


el eRe 
.e 
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Assuming that it is desired to obtain the expression for the primary | 
error in 22, the desired expression, in accordance with Equation (11b), 
is as follows: 


Cie Co C2 
da. = é + e2 + es (15) 
x = 0 
Aoi a3 an 13 Qi a3 
where (OFF =>— , Cx = ; and Cre = 
Q31 33 Q31 33 Qa Ae3 


Consider next the following unit system, formed in the manner 
specified in the foregoing: 
Qn%1 + Ant, + da%3 = 0 
Q12%1 + Ao2% + ars = —1 F 
341 + Ae3%2 + 3323 = 0 : 


A solution of these equations gives the following values for the 
unknowns: 
1 


ty = — 


D 


=i 


D 


fos 1 
;and 7; = —— 


D 


Goi 31 Qi 31 Qi =A 


) 


Q23 33 Qizs 33 a3 Ae3 


By referring to Equation (15) it is evident that the foregoing 
values for % .... %; are identical, respectively, with the coefficient 
of e,.... 3 in the expression for 82». . : 

This method of obtaining the expression for the primary errors 
in the unknowns has considerable practical value, inasmuch as it is 
applicable not only to all ordinary systems of equations, but also toa 
class of equations of great importance in engineering problems, 
namely, those which are solvable by successive approximations. In 
the former case, the corresponding unit systems may be solved by 
determinants or by elimination ; in the latter case, by iteration or 
by converging increments.* 

Another important class of equations, which may be regarded as 
a special case coming under the present method, are those possessing 
the property of diagonal symmetry. The Clapeyron equations in 
mechanics and structural analysis (including the equations of slope 
and deflection) and the normal equations in Least Squares are 
generally of this type (see Example No. 11). 

Referring to Equations (1) to (1n), inclusive, Section 5, if the 
initial values of the coefficients are such that the coefficient in the 
i-th row and the k-th column is the same as that in the k-th row and 


“Solution of Equations in Structural Anal sis by Converging Increments,” by G. H. Dell 
Proc. A.8.C.E., October, 1938. or ging e y ell, 


a 
: 
, 

, 
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the i-th column for all values of 7 and k from 1 to n, a system of 
diagonally symmetrical equations is obtained.* In determining the 
primary errors in the unknowns of equations of this type, it is ob- 
viously unnecessary to make any re-arrangement of the coefficients 
when forming the corresponding unit systems. 


Example No. 3. Derive the expression for the.primary error in 2°, 
Example No. 1, by use of the unit system. 


Solution: The unit system is as follows: 
Spchey ie at 1 
3 — Sy +22= 0 
ape op e708 
A solution of these equations by determinants or by elimination gives 
the following values for the unknowns: 


-1_—-9 4) 


x= —j y = ——; andz = —— 
32 32 32 


Hence the expression for 6z is 


(é1 + 9e2 + 21es) 


r= 
— 32 


11. Finding the Maximum Primary Error.—In order to determine 
the maximum primary error 6%m in a given unknown zm, after the 
expression for 6zm has been obtained in terms of @.... én by one of 
the preceding methods, this expression is expanded in accordance 
with Equations (7) and (4), Section 7. There is thus obtained an 
expression in which the errors Aki .... Ak, occur linearly. In the 
case of interdependent coefficients, a re-arrangement of terms is 
necessary in order to combine the multipliers of a given Ak into a 
single error-coefficient. 

In order to obtain the maximum positive value of the primary 
error in 22, the signs of the increments Ak; .... Ak», should be the 
same as those of the corresponding error-coefficients, and, vice versa, 
in order to obtain the maximum negative value of the primary error 
in x°, the signs of the increments Ak; .... Ak» should be opposite 
to those of the corresponding error-coefficients. After the signs of 


*Note: While, in order to be diagonally symmetrical, the relation ak = aji must necessarily hold, 
the coefficients aj. and az; (considered as variables) are not necessarily interdependent and may become 


unequal. 


lias 


Fb 7 S fs “> he : , oe an - 
= als z 7 ae Re Ve oe “—- r tv - et ee 
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the Ak’s have been thus determined, by inspection, the num 
“EG of the maximum primary error is readily calculated. 


_ Example No. 4. Assuming that the coefficients in the equations 
of Example No. 1 are subject to errors of +10 per cent, determine 
the maximum primary error in 2°. 

Solution: By the methods of Sections 8, 9, and 10, the expressio n 
for the primary error in x° was found to be 


i : a 


where é&: = x Aky + y° Aki. + 2°Aki3 
é€2 = 2°Aka, + y? Ake 7 2°Ako3 ~ 
and es = x°Aksi + y? Aks2 a 2°Akss 


Upon substituting the initial values of the unknowns, the expression — 
yy for the primary error in z° assumes the following expanded form: 

2Aku + 3Aki2 — 0.05Akis % 

on Vang + 9(2Ako; + 3Akee ~— 0.05Ake3) “aq 

+ 21(2Aks: + 3Ak32 — 0.05Akss) ‘ 


Fy 


= 
s 


- In order to make each term in this expression a positive quantity, — 

7 the following signs must be used with the various Ak’s: > 
ey = ae = 3 
=> = aie 
= a id 


and upon giving the Ak’s the specified magnitudes of 10 per cent, the 


maximum positive value of 62 is 
f 


2(—0.2) + 3(—0.3) — 0.05(0.3) 
da (max.) = a5: + 9[2(—0.1) + 3(—0.5) — 0.05(0.2)] 
+ 21[2(—0.1) + 3(—0.2) — 0.05(0.1)] 


1 
aie rcs [—1.315 + 9(—1.710) + 21(—0.805)] 


= +1.050. 


. > 


RS IN i SIMUL T ANE oe) ta eos vA: Pi 7 19 ’ 


a” p- 


re ie. : , 
For a maximum negative primary error, the signs of the Ak’s 
d all be reversed, as follows: 


+ + = 
+ + - 
+ + - 


and the maximum negative value of 5z is — 1.050. 


a | 


12. Primary Errors in the Unknowns of Non-linear Equations.— 
1 he following functions will be used to represent a system of n 
linearly independent non-linear equations in n unknowns: 


Py — Orn +.... + baXin +-.-. + bunXn —m =0 (17) 


os 


Pg Unt Amt It -s oot OneAen 2s s). 1 CaaXnn — Vm = 0 (17m) 


Fy, = bu Xu +....+ OnmXam +.--- + DanXnn — Vn = 0 (17n) 


In a typical equation, namely, 


FF; = baXat+.... + DimXim +. ... + dinXin — 15 = 0 


the quantities Xi; .... Xi, are functions of one or more of the 
unknowns 2 .... 2,; and the coefficients bi .... bin are functions 
of one or more of the p arbitrary constants ki .... k, occurring in 


the given system, which are subject to the errors or variations 
Ak, .... Aky, respectively. 


The quantities Xi; .... Xin may be represented by the typical 
function X;;(z1, ...., %n), whence 
OX ij OX 5; 
5X i; = aX 3; — ( 621 a eealale + ire) (18) 
° Kol In ° 


A typical coefficient may be represented by the function };;(ki, 


. , kp), whence 
(<2 oa he 00;; ee ) 
a aki 1 oe ee dp Dp 


(19) 


ie) 


6b:; = db;; 


The same assumptions which were made in Section 7 regarding 

the magnitudes of the errors in the arbitrary constants will again be 

used. It will also be assumed that the original equations have been 
solved for the initial values z?.... 2n. 


: 


a : 


_ hee ee A, 


| ae 2 —— = a ph AP 
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From ete ae ClTY tor Cl 7a); inolusive. the total differen 
GPa. Sof OF aye. . a5 GF, are as follows: 


aF aF 
Tp eNotes ey ise eee 
O21 0 


m 


Te Xpdbyot ON abe pee eee 


OF n OF. 
3 dim te... et 


Lm Ons 7 


OF nm 
ha dxa+....t+ 
Ox 


m = 
din 


+ Xmdbm +....+ Xmmdbmm +....+ Xmndbmn = 0 (20m) | 


— 


OF , oF ,, 2 WES 
= da+....+- dim t+....4+ dXn 
0x1 0 


Lin OLn 


eX midOni. bv sn st Xam @Dam na be (20n) 


From these equations the corresponding relations between the : 


errors in the coefficients and the primary errors in the unknowns ~ 


: 
may be obtained by substitution of the initial values, bf, 22a Von . 
and z?.... 2°, as follows: s 
aFy oF AF, n . . 
duit... .+ Otm+....+ 62n+ > X?,5b,;,=0 Cam 
O21 |o Lm lo Oita 1S q7=1 a | 
3 
files dai+ neces d2m+ Tee 3 X° 8b, =0 (21 ‘ 
ieee Unaiaerete a mjo0m; = Pe 
ork 02m \o OXn ‘ pd / Ce g 
¥ 
& 
oF, OF, OF, 
duit... .4+: bim+... + ito + DXe, 6b,,;=0 ‘(inks EB 
Ox oO Osta i) Xn \o 


As in the case of linear equations, Section 7, it will be desirable : 


for the present to employ an abbreviated evmbal for the terms — 


indicated by the summation sign in Equations (21) to (21n). For 
this purpose, let 


7 be X?,6b;; (22) | 


The following linear equations may then be written in place of — 


Equations (21): 


@iy 


(21m)’ 


; OF, 
—|6m+....+ btm t+....+ btn = —E, (21n)’ 
OX, Crs Ove ° . OXn o 
From these equations, the expressions for the primary errors 621... . 


5x, may be obtained by the methods discussed in Sections 8, 9, and 10. 
It should be noted, however, that the coefficients to be used are those 
indicated in Equations (21), and (21)’, and not those in the given 
non-linear equations. 

In order to derive working expressions to be used in connection 
with the method of determinants, it will be noted that the deter- 
minant of the coefficients of ape unknowns 62,.... 62, in Equations 
(21)’ is the Jacobian 


OF; OF, OF, 
an, Kg ee rs aa, E mise 6 ax, i 
ore OF n OF n 
J = See Sn Ae (23) 
Ox ° OXm (s) Oly ° 
OF, OF, OF, 
ae . ep te Bice i cis ¢ Se : 


The primary error in a given unknown, 2%, is then given by the 
following solution: 


OF, OF, 
0X1 To) ohien °o 
1 OF m OB 
62m = —— aie (24a) 
—J 0X1 |o OtneG 
oF, OF , 
02, oO OXn ° 
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Ye 


" et! a , 
, & Dee bes wd 
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ae 
4 2 hk 

Denoting the co-factor of the element in the i-th row and the 
m-th column of the Jacobian by the symbol Gin, and developing the 
foregoing determinant with respect to the elements of the m-th 
column, the primary error 6z, may be expressed in the following 
abbreviated form: 


ip se : 
82m = —— >> GimE; (24b) 


— t=1 


: The routine to be followed in practical problems concerned with _ 

i the calculation of the primary errors in the unknowns of non-linear | 

if equations is quite similar to that described in connection with linear | 
equations. After the coefficients of H,....H, in the expression for 
5am, Equation (24b), are obtained, the expression is further expanded ~ 
in accordance with Equations (19) and (22). This procedure results 
in an expression in which the errors Aki .... Ak, occur linearly, and 
it is thus possible to determine by inspection the proper signs which _ 
should accompany the Ak’s in order to obtain the maximum primary 
error in a given unknown. 


Example No. 5. The equations of two intersecting curves are as — 
follows: . 


Fy = ky? & box do eg = Oe 6 (25) | 
F. = (kax)? + kaksxy + (ksy)? — 102.25 = 0- (26) a 
With the following initial values 
ki = 1, kp = 2, kf = —8, kf = 1.5, and ke = 2, Se 
the original equations are: 
(x°)? + 2° — 3y° = 29 
2.25(x°)? + 3x°y? + 4(y°)? = 102.25 ; 


and the curves intersect in four points, one of which has the values _ 
x = 5, y° = 2. ; 
It is required to determine the maximum primary error in y° at — 


this point due to possible errors or variations of 10 per cent in the 
arbitrary constants k,.... ks. 


iS | 2kex°e + ke ke 
2(kg)Pae + kgkgy? §—-2(ks)?y? + kghkgae 


[iz -3 
| es rep ey = 457.5. 
_ By Equation (24a), 
_ aFy 
| BE, 
kal ax |. -- a ; 
—J | oF: —J |28.5 E, | 
Ox |o : 
et 28.5EH, — 12E, (27) 
457.5 


= © 


The values of HZ; and E:, according to Equations (19) and (22), ; 
are the differential changes in F; and F», respectively, due to varia- 
tions in the arbitrary constants ki .... ks, namely, . 


7 Ey => (z°)?Aki = r°Ake +y°Ak; = 25Ak, + 5Ake — 2Aks3 
Ey =[2k2(x°)? +koxey?|Aks+[2k3(y°)?+hixey?|Aks = 95Ak.+31Ak; 


Referring now to Equation (27), it is seen that, in order to obtain 
maximum positive Sy, the errors in the arbitrary constants should 
_ have signs as follows: 


en iniak ss 0:2) AR = 40.3 Akg = 0:18 Ake = ~0.2 


_ Upon substituting these values in the foregoing expressions for E; 
and E> it is found that EZ, = 4.1 and E#, = —20.45; and, from 
Equation (27), 
28.5(4.1) + 12(20.45) 
457.5 


= +0.792 


dy (max.) = 
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The calculation of the corresponding primary error in 2° is as 8 


follows: 2 
OF; 
Nea es ay 1). (2; . 
ie = ——— —— % 
—J OF, —J Ey, 31 ‘ 
E, 
dy 
lz 3H 31(4.1) — 3(20.45 
—457.5 — 457.5 


It is thus seen that the coérdinates of the point of intersection of 
the two curves are changed from (5, 2) to (4.856, 2.792) by the 
assumed errors in the arbitrary constants. Since only the primary 
errors were considered, the foregoing results should be regarded as 
approximations. 


13. Effect of Errors in Right-hand Terms. (a) Linear Equations — 
If, in addition to the errors in the coefficients of Equations (1) to 
(1n), inclusive, the w terms are subject to the errors or variations 
Api... . Apn, the total differentials in Equations (5) to (5n), inclusive, 
will include the terms —du, .... —dy,, and the following relations 
will be obtained in place of Equations (6) to (6n), inclusive: 


At 6%, +... + OF, 5%_ +... + 08,52, + >) 2A; — Aum =0 (28) 
fest 


G50, +... + On dm +... + 09,,52,+ > 2°Admj—Aum=0 (28m) | 
j=1 


An0%, HF -- + Fim dm +... + 2,80,+ > 29AGnj—Atin=0 . (28n) | 
j=l f | 


The expression for the primary error in rz» then becomes 


loners 
0Lm = —— DS Cim(es — Ap) (29a) 
—-D ia 


If the right-hand terms alone are subject to errors or variations, 
the e terms in the foregoing expression are zero. In this case, 672m, 
» 6"%m in Equation (3b), Section 6, are equal to zero, and the 
sere error and the true error are tential, namely, 


pe 


_ 


oo 


on ak ee ees 7 , : ; 
_ ERRORS IN SIMULTANEOUS EQUATIONS | 25 


es 
Ax$, = 62m = ay >> CimApi (29b) 
i=1 


(b) Non-linear Equations. If, in addition to the errors in the 
coefficients of Equations (17) to (17n), inclusive, the » terms are 
subject to the errors or variations Av, .... Av,, the total differentials 
in Equations (20) to (20n), inclusive, will include the terms —d» 
+++. —dv,. The vy terms may consist of one or more arbitrary con- 
‘stants occurring linearly or non-linearly. In any case, the true maxi- 
mum errors Ay, .... Av, may be calculated. 

_ The following relations are then obtained in place of Equations 
(21) to (21n), inclusive, 


OF, OF, OF; as 
dti+....+ bim+....+ 52,+ >> X2,5b,; — An = 0 (30) 
021 |o Zm |\o Xn \o j=1 
OF , OF, OF n n 
6ai+....+}—| 6tm+....+—| 62,+ >> X°2,8b,,; —Avm =0(30m) 
Z1 Io Olmlo OL; lo ph 
aF,, aF , aF 0 
dai+....-+ Olg tie atk 6tn+ >) X°6b,;—Avn=0 (30n) 
1 |o OLmia OLa\» jZ1 


The expression for the primary error in z, is then as follows: 


7 eae 
i it. = —— a Gim( Ei; — Avi) (31a) 
= fi Pm 

; If the v terms alone are subject to errors or variations, the # 
terms in Equation (31a) will be zero, and the expression for the 
primary error becomes 


{ 2 
62m = ire oe GimAvi (31b) 


i=1 


An important use of the relations expressed by Equation (31b) 
occurs in connection with the solution of non-linear equations, as 
_ will be explained in Section 16. 


Example No. 6. In connection with the equations of Example 
No. 5, find the maximum primary error in y° and the corresponding 
primary error in 2°, due to errors of 10 per cent in the right-hand 


; terms. 
P 


: 
f 


i 


t ib), a er a iy 


; i oa : Ai a Pe om 33 | 
28.5, 28.50 + 12An we 
Wea oe Oiea re 457.5 "0 
ee 
See Av; 33 
ie eek Ui, Ary 
oe TNA GT <2) lah hye LOTR 
ee ‘a 


In order to obtain the maximum primary error in y°, Av and An 
should have opposite signs. Hence : wal 


a 28.5 (2.9) + 12 (10.225) 
2 457.5 
| 31 (2.9) — 3 (10.225) 
j oS 6x ae =e Se ee See 
: 457.5 


lI 


dy (max.) = +0.45, and 


= +0.13. ‘ 


14. Partial Differentials and Partial Derivatives.—It has been seen 
Fe that the expression for the primary error of a given unknown, in- 
the case of both linear and non-linear equations, can be put in the 

following form: 
_ btm = cr Aky +.... + ch Ak, (32) 


~ 


“or, using the differential notation, 


dim = 01 diy +o. a0 + Cy Ob, (33) 


where dz,, is defined as follows: 


diy = dig th fc. ee ae 1 (34) © 


In Conon (32), the coefficients of Ak, .... Ak», which will be 
designated as ‘‘error-coefficients,”’ are hinted ‘ the initial values” 
1 GER Steg OMe Ns Ns Bg eek 

The relation exorened by Equation (33) is a more general one, 
in which the error-coefficients c, .... cp are functions of any valve q 
of ky.... kp and a .... 2, which are satisfied by Equations (1) to. 
(1n), clase or by Equations (17) to (17n), inclusive. A com- 
parison of Raunone (33) and (34) shows that the terms c, dk, 


~ dh = —— dky, respectively; 
Pp 


bso ) that the error-coefficients c; .... cp are the partial derivatives 
: a respectively. 4 
_ In Section 17 it will be shown that the determination of the true 
maximum error involves the investigation of the signs of these partial 
derivatives for various values of the arbitrary constants (and for the 
corresponding values of the unknowns). 

In regard to the effects of errors or variations in an individual 
arbitrary constant, k;, of a system of simultaneous equations, the 


primary error in z, is given by 


Ak;. Using initial values, 


Equation (34) may be written as follows: 


OD 
Aki +....+-—— 


° p 


Sim = Ak, (35) 


oO 


Bi: 


from which it is seen, upon comparison with Equation (32), that the 


primary error Ak; caused by a variation in an individual arbi- 


°o 


_ trary constant may be obtained from the general expressions which 
_ were derived in Sections 8 and 12 by placing the remaining Ak’s 
equal to zero. 


15. Probable Errors.—In the preceding section, considering the 
; unknowns 2 .... 2, of a system of simultaneous equations as func- 
tions of the arbitrary constants, it was pointed out that the primary 


error 62m is equivalent to the sum of the primary errors Akiecer 
} 1 lo 

OXm : : : 

Ak,, obtained by allowing k,...., to vary one at a time. 

0 
- p\o 
In a similar manner, if the arbitrary constants ki .... kp are 
subject to the probable errors +Ak; .... +Aky, the corresponding 
i 

probable errors in a given unknown, z;,, are + FN gel es 
; Okan : 1 Jo ‘ 
B+ Aky, respectively. 
m Ok, |o 


The probable error in x°,, due to these combined effects, is found 
by taking the square root of the sum of the squares of the individual 
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effects;* or, denoting the probable error in 2;, by the symbol )z*%,, z 


OLm 2. OXLm 
Ag ahs) +... +( 
Ok; 


ak, 


It will be noted that the application of Equation (36) is facilitated 
by first writing the expression for the primary error, as given by 
Equations (11b) and (24b), and then reducing to the form indicated 
by Equation (32). 

If the effects of probable errors in the right-hand terms are to be 
included in the calculation of the probable error in z;,, Equation (36) 
should be extended by including under the radical the sum of the 
squares of the corresponding terms occurring in the expression for 
52m according to Section 13 (see Example No. 8). 


oO 


Ap) (36) 


Example No.7. Let it be required to calculate the probable error 
in x° of Example No. 4, assuming the arbitrary constants to contain 
probable errors amounting to +5 per cent. 


Solution: In Example No. 4, the expression for 6x was found to 
be as follows: 
(2Aky, + 3Akiy2 — 0.05Ak3) 


62 = ar + 9(2Aka: + 3Ake2, — 0.05Ak3) 
+ 21(2Ak31 + 3Ak32 — 0.05Ak33) 


In accordance with Equation (36), the probable error in x° is then 


1 agri e (0.2)? + (0.45)? + (0.0075)? 2 
~ 39 | + 81[(0.1)? + (0.75)? + (0.005)2] = +0.298 ~ 
+ 441[(0.1)? + (0.3)? + (0.0025)?] 


Example No. 8. Find the probable error in y°® of the equations 
in Examples No. 5 and No. 6, assuming that the arbitrary constants 


and also the right-hand terms are subject to probable errors amounting 
to +5 per cent. 


Solution: In Example No. 5, the primary error in y° due to errors 
in the arbitrary constants of the coefficients was as follows: 


ae 1 —28.5(25Aky + 5Ak, + 2Aks) 
—457.5 
+12(95Aks + 31Aks) 


*Art. 3, Chap. III, Geodesy, by A. R. Clarke. 


oh eee! ene aaa 


SS ES ee eee 


ay, 


shown in Example No. 6, is 


An and Av; and proceeding in accordance with Equation (36), the 
following computation gives the required probable error in y°: 


xye = 2 _ /812.25{(1-25)? + 0.5)" + (0.3))] 
i | 45755 + 144[(7.125)? + (3.1)?] = +0.274 


+ 812.25(1:45)? + 144(5.1125)? 


IV. True Errors 


16. Methods Available for Calculating the True Error—The true 
error in the unknowns of a system of simultaneous equations, due to 
errors of definite magnitude and sign in the arbitrary constants, can 
be calculated in three ways, as follows: 

(1) By writing and solving the revised equations. 

(2) By using approximate values for the unknowns in the revised 
equations and then varying the right-hand terms and calculating the 
corresponding corrections to the unknowns (as in Newton’s Method). 

(3) By Taylor’s series, involving the calculation of the errors of 


_ higher orders. 


These methods are all of value in investigating the true maximum 
error. The first is best adapted to linear equations, and the second 
to non-linear equations. The third method is more limited in its 
usefulness as a tool for numerical calculation, but has considerable 
analytical value. A brief explanation of the various methods follows. 

First Method. Upon substitution of the modified coefficients, 
Qi, .... Qn, in Equations (1 to (1n), inclusive, the revised equations 


are 
ane +... tant +... + aint, — wi = 0 (37) 


Git ee + Onmtm eee Omnntn — p= 0 (37m) 


Gee eel Cynbw «oe Hh Oantn — bp = 0 (37n) 


The primary error in y° due to errors in the right-hand terms, as 


i eal a 


1 
dy = ——— (28.5Ap; — 
7] —A575 ( Vi 12Ap2) 
Upon substituting the assumed values of Aki .... Aks and of 


Rs 
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The true error in a given unknown, zm, is then given by Atm = Im — 
z®,, where 2, is obtained by solving the revised equations, and 2m 
is the initial value of the unknown. | 

Second Method. In considering the application of this method 
to linear equations, let it be supposed that the unknowns zi .... @n 
in Equations (37) to (37n) are approximately given by a1 .... 2n, 
respectively. It is generally convenient to make use of the primary 
errors in obtaining such approximate values, in which case 1 = 2} 


+ 621, v2 = 28 + 6x2, ete. Next, suppose that values of ui.... pee 
are computed from the relations, 
Ait +e... + Oimtm +... Ht Ointn = MI (38) 
Odi +... + Onminm +. ... + Onntn = bm (38m) 
nit fee HOmTm tos. bt Oneten = oh (38n) 


If the right-hand terms of Equations (38) to (88n) are now given the 
variations Au. = pf — pi, Aus = uo — ps, etc., the method described 
in Section 13 may be used to calculate the corresponding variations 
A’x,.... A’, in the approximate values of the unknowns, 21... . 2n} 
and, as stated in Section 13, these variations will be the true errors, 
or true corrections, when dealing with linear equations. The true 
error in zy, is then Az’, = 62m + Atm. 

The application of this method to non-linear equations is similar 
to that explained in the foregoing, with the exception that when the 
right-hand terms »;.... , are given the variations Avy, = vf — v, 
Av, = vz — v2, etc., the corresponding variations obtained in accord- 
ance with Equation (31b), Section 13, are not the true errors in the 
unknowns; consequently, the process is repeated until the corrections 
to the unknowns become negligible. 

It will be noted that this procedure affords a convenient method 
of solving non-linear equations after approximate values of the un- 
knowns have been determined. 

Example No. 10 illustrates the use of this method in connection 
with non-linear equations. 

Third Method. If the expression for the primary error of a given 
unknown, x, of a system of linear equations is changed in such a 
way as to replace the unknowns zi .... x2 by the primary errors 
6a, ... .6%,, respectively, the resulting expression represents the value 
of the secondary error, 6x». . Likewise, if in this expression for the 


secondary error, the quantities 62; .... 6x, are replaced by 2. 

x, respectively, the resulting exniesion represents the value of the 
third order error, 6*z,. In a similar way, one may obtain each term 
in the Taylor’s series, Equation (3b), Section 6, after the errors of 
the next lower order have been calculated. The method of deriving 
these results is given in the Appendix. 

This method of calculating true errors appears to be chiefly useful 
in studying the effect of variations in individual arbitrary constants, 
or in small groups thereof, in connection with the determination of 
the true maximum error. 

The expressions for the errors of higher order in the unknowns of 
non-linear equations become increasingly complicated, and are there- 
fore of insufficient practical value to warrant their inclusion in the 
present treatment. 


17. The True Maximum Error.—In investigating the true maxi- 
_ mum error in a given unknown of a system of simultaneous equations 
when the assumed errors or variations Ak, .... Ak, in the arbitrary 
constants ki .... k,, respectively, are assigned definite limits, the 
usual concern is the determination of the correct signs to accompany 
the Ak’s. In practical engineering problems, the true maximum 
error requires, as a general rule, the use of the assumed maximum 
value for each of the Ak’s; in other words, a true extremum, as 
treated by the methods of maxima and minima, will seldom be a 
condition affecting the true maximum error. Consequently, the 
7 present discussion will deal principally with cases which are sus- 
__eeptible to treatment by elementary methods. 

7 Since a given unknown zm of a system of linearly independent 
equations is fully determined for definite values of the arbitrary con- 
stants, it will be convenient to consider z,, as a function of the inde- 
pendent variables, ki .... k,, as follows: 


Tm = J(hi, .~..-, Kp) (39) 


_ The true maximum error will be obtained by making |zm — zn| as 
great as possible without permitting the Ak’s to exceed their assumed 
maximum values. 

Let Zn — 2°, denote the maximum positive error in tm, and let 
2m — x°, denote the maximum negative error. Then, by means of 
Pe oropriate variations in the magnitudes of the Ak’s, x» may be made 
to vary continuously throughout the region 2m < Lm < Xm. As the 
ensuing arguments will apply equally well to the maximum positive 


32 ILLINOIS ENGINEERING EXPERIMENT STATION 


error and to the maximum negative error,.it will be sufficient to 3 
confine the discussion to the region x7, = %m = Um. ‘ 

The procedure that will be used in arriving at the true maximum ~ 
error may be outlined as follows: Using the assumed maximum 
values of the Ak’s in conjunction with such signs as appear most 
likely to lead to the maximum error in 2%, the true values of the 


unknowns, v,.... 2, are calculated. In this way a point 2, in the 

region under consideration is obtained, and the signs of the partial 
Of OLm ; : : 

derivatives ae are then investigated. If these signs are 
Oky dk» 


such that a further increase in x, — x°, can be obtained by increasing 

each of the Ak’s beyond their assumed maximum values, the dif- 

ference zx,, — x°, is the true maximum error. If, on the contrary, 
OXLm OLm 


aie eae 
point are such that a decrease in the absolute value of the corre- 
sponding Ak’s will produce a further increase in z,, — x?,, the signs 


at this 


the signs of some of the partial derivatives, say 


of Ak, .... Ak, are reversed, and a new set of modified unknowns, 
z, .... 2), are calculated. The signs of all of the partial derivatives 
OXm OXm ee. ‘ : ae 

fy. are again investigated; and by applying the criteria 
Oki Ok» 


previously described it is learned whether or not the difference 
x’ — 2° is a maximum. If the condition of an extremum is not in- 
volved, the true maximum error will be obtainable by means of the 
simple process described. 

If an extremum should be involved, or, in other words, if the 
maximum error occurs when for some of the arbitrary constants the 


Xm 


partial derivatives have zero values, it becomes necessary, in 


general, to use values for the corresponding Ak’s which are less than 
the assumed maximum values. In such cases, it is still necessary 
to use the assumed maximum values of the Ak’s for those arbitrary 


O2m 
constants for which Re does not vanish. 


The process of determining the signs of the partial derivatives 
Otn OL 5 
oF at the point 2, is most readily performed by taking 


Oki Ok, 
the expression for 6%m (the primary error), and substituting therein 
the modified values of the arbitrary constants (k,.... k’,) and of the 


/ / s . aes, 
unknowns (2, .... x,) in place of the corresponding initial values. 


ERRORS IN SIMULTANEOUS EQUATIONS oo 


pach an expression will, for convenience, be referred to as the 

“criterion expression,” or “the expression for 6z/,.””. The error- 

coefficients in the criterion expression will be the values of the partial 
Olm OXm 


Ok, dk 

A brief examination will now be made into the reasoning under- 
lying the procedure which has been outlined for finding the true 
maximum error in xm. Referring to the function indicated by Equa- 
tion (39), it will be assumed that x°, has been changed to another 
point, z,, in the region x}, = tm = Xm, by means of certain incre- 
ments Ak; .... Ak, which have been given such signs as are required 
for a maximum positive primary error, 52». If the values of the Ak’s 
are taken sufficiently small, the signs of the partial derivatives 


Len Gin 1 o. , 
ee. will remain unchanged; accordingly, in order to 
Ok, 0k, 


secure a further increase in the difference x, — 2°, that is, in the 
true error, Az,,, it is sufficient to increase the value of each of the 
increments Ak; .... Ak, by a small amount. Furthermore, in the 


derivatives at the point in question. 


= remains unchanged 
in sign throughout the interval 2), = tn = %m, 1t is necessary to use 
the assumed maximum value of Ak; in order to make the difference 
Lm — Lm as large as possible. 

Suppose now, following the procedure outlined earlier, that the 
unknowns x} .... 2° have been changed to 21 .... x;, respectively, 
by using the assumed maximum values of the Ak’s in conjunction 
with the signs required for making the primary error, 6z,, a maximum. 
From the preceding discussion it is seen that if none of the partial 

OG, OLm 
derivatives ——.... 
1 P 
varied from z°, to z},, the maximum error in m, iS tm — Lm = Lm — Lm. 

Let us now consider the case where some (but not all) of the 

On OXLm 


Ok», Oky 
changed from 2°, to 2. Up to the present time, those Ak’s for which 
the partial derivatives have undergone no change in sign have ren- 
dered their full contribution towards obtaining a maximum error in 7», 
and it will be advantageous, therefore, to treat these arbitrary con- 
stants as ordinary fixed constants while investigating further changes 
in the remaining arbitrary constants, and to consider the function 


Pepe GHA wo cs gh) 


case of every arbitrary constant for which 


have undergone a change of sign as 2» 


partial derivatives, say, , have changed sign while 7» 


in x, the value of am — 2m will be further increased if the increments 


increase. Beyond this stage, by reversing the signs of the Ak’s in 
question, and assigning them sufficiently small values, it will be — 
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in the region 7, = %m = Em. Since the signs of the partial derivatives 
Oty Om 


SaaS 
those of the error-coefficients in the expression for the primary error | 


at the point z,, are now assumed to be opposite tom 


Ak, .... Akg, each or all, are reduced by a sufficiently small ane | 
from Picts assumed maximum values. Furthermore, if these partial 
derivatives undergo no further change in sign while Ak, .... Ak, are — 
gradually reduced to zero, the value of zt, — 2m will continue to 


possible to cause the difference x, — x to increase still further; and — 
if, as the Ak’s approach their assumed maximum values (their cn 
DLin ¥ 


having been reversed) the signs of the partial derivatives ak, sa 
On . 


ak do not again change, the maximum error will be obtained by 
k 4 


using the assumed maximum values of Ak, .... Aky. 

It is thus seen that if, following the recommended procedure, a — 
new set of modified unknowns, zj/ .... 21, is calculated after re- — 
versing the signs of Ak, .... Ak, (whose magnitudes are taken equal 
to their assumed maximum values), and the partial derivatives 
OD OXLm 
Ok, a Ok» 


increase in %», — x» Will be obtained by increasing each of the Ak’s 
beyond their assumed maximum values, the true maximum error is 


are all of such sign, at the point 2)’, that a further — 


given by rm — tm = Im — Lm. ~ 
It may be found necessary in some instances to repeat the above ~ 
process for a third set of unknowns zj’’ .... x’, before the true 


maximum error can be found. ‘ 


In some problems it is possible to foresee a change in sign of some of — 


; ' : OLm OD 
the partial derivatives : 


before the true values of 
1 p 
v,.... 2, are calculated. For example, if, when signs are used with 
the Ak’s as required to obtain a maximum primary error in 2°, . 
approximate values of z1.... 2), are found from a computation of — 
the primary errors 62, .... 6x,, and if it is thereby indicated that 
some of the unknowns may change sign, it is seen, by referring to — 
Equation (11b), that the signs of some of the quantities «’Aa;; in 
the corresponding criterion expression may likewise change. It is 
therefore advisable, before trying for a true maximum error, to calcu- 


r vat » ’ Ve : o 
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> the primary errors in all of the unknowns, and, if it appears 
bable that some of the error-coefficients will change sign, to re- 

: se the signs of the appropriate Ak’s before computing the apparent 
m aximum error. 

Finally, in regard to cases involving an extremum, it may be © 
said that if, after the difference zx, — x, has been made as large as 


possible by using the assumed maximum values of the Ak’s with those 


Lm 


\ % . . *. ? . . 
arbitrary constants for which the partial derivatives do not 


vanish, there remain some arbitrary constants for which the partial 


OL 


derivatives 


do vanish when the absolute values of the corre- 
sponding Ak’s are less than their assumed maximum values, the 
methods of maxima and minima must ordinarily be used to deter- 
mine the exact values which those Ak’s require in order to secure 
the true maximum error. It is difficult to formulate definite rules 
setting forth a general procedure for dealing with such cases without 
expanding the present study far beyond the limits contemplated. 
_ However, it will be observed, first, that the procedure previously 
_ described enables one to determine certain regions in the neighbor- 
hood of %, in which the possible existence of an extremum will be 


indicated by a change in sign of site when the sign of the accom- 
panying Ak is reversed; and, second, to determine for which of the 
arbitrary constants the assumed maximum values of the Ak’s should 
be used. 

Example No. 9. From the equations used in Examples Nos. 1 
and 4, let it be required to find the true maximum positive and 
- negative errors in x° due to the assumed maximum errors of 10 per 
- cent in the arbitrary constants. 

Since there are nine arbitrary constants in the given system, the 

total number of ways in which the signs of the specified errors can 

- be combined is equal to 2° = 512. The problem involves finding 

two particular combinations, namely, those required for producing 

the maximum positive error and the maximum negative error, 
respectively, in 2°. 

The signs to be taken in computing the apparent maximum error 

in x° are those which are required to render the primary error a 

maximum, unless it can be foreseen that some of the error-coefficients 


Ox Ox 
Oky ” Okie 


, ete., undergo a change in sign while z is changing from 


—— 
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z° to x’. This will be the case if the sign of 2’, y’, or z’ is different — 
from that of x°, y°, or 2°, respectively. In this event, for purposes — 
of computing the apparent. maximum error, the signs of those Ak’s 
whose coefficients contain the unknown which has changed from posi- 
tive to negative, or vice versa, should be taken opposite to those . 
required by them for making the primary error in z° a maximum. 
Therefore, in order to make a more intelligent estimate of the signs 
required for the true maximum error, it is expedient to first calculate 
the primary errors in all of the unknowns, as in this way it is generally 
possible to determine whether a change in sign occurs in any of the 
unknowns. 

The expressions for 62, dy, and 6z were found to be as follows: 


il 
oS 30 ( e1 + 9ey a 21e3) 


1 
éby = —— (3e, — 5 —-— ee 
y sae 1 2 3) 


1 
6z = aan! (Je, a RCT a 13¢e3) 


where, for maximum positive 6x the values of ¢1, e2, and es were found 
to be —1.815, —1.710, and —0.805, respectively. Hence, when 
6x (max.) = +1.050 


1 
by = rare (—3.945 + 8.550 + 0.805) = —0.169 


1 
and éz = = (—9.205 + 1.710 + 10.465) = —0.093 


The modified values of the unknowns will be approximately as 
follows: 


x’ = 2.000 + 1.050 = 3.050 
y’ = 3.000 — 0.169 = 2.831 


z’ = —0.050 — 0.093 = —0.143 


therefore, since none of the unknowns have changed sign, the apparent 
maximum error will be calculated by using with the Ak’s the same 
set of signs as was used in computing dz, namely, 
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The revised equations are then 


1.82’ + 2.7y' + 3.32’ = 12.85 
0.92’ — 5.5y’ + 2.22’ = —13.10 
0.92’ + 1.8y’ — 0.92’ = 8.05 
for which the solution is as follows: 
D'= 1.8(0.99)+ 0.9 ( 8.37)+0.9 ( 24.09) = 30.996 
a : [12.85(0.99) —13.10( 8.37)+8.05( 24.09)]= 3.129 


~ 30.996 


: ; 
y’ = 309061 12:85(2.79) — 13.10(—4.59)+8.05(— 0.99)]= 2.839 


z= 5 [12-85(6.57) —13.10(—0.81) +8.05(—12.33)] = —0.1362 


The apparent maximum error (positive) in z° is therefore 3.129 —2.00 
= +1.129. 

Test: As previously stated, in order to determine whether or not 
the foregoing result is the true maximum (positive) error in 2°, it is 
necessary to see whether or not the products, or differentials, 

/ / 
aie Aku, ae Aki, ete., are all of like sign. This process is equiva- 
Oki Oki. 
lent to taking the revised equations and writing out the expression 
for the primary error in 2’, and investigating the resulting signs of 
the separate terms, using for the Ak’s the same signs as were em- 
ployed in the revised equations from which the apparent maximum 
error was computed. As previously observed, the expression for 
the primary error in 2’ is obtained from the determinant solution 
of x’ by substituting ei, ¢, and e3 in place of the uw values, 12.85, 
—13.10 and 8.05, respectively, and dividing by —1; hence, 


1 / 
dx’ = ————— (0.99, + 8.37¢, + 24.09e;) 
— 30.996 
where e, = 3.129 Aky, + 2.839 Aki, — 0.1862 Akis 


e3 = Bh WA Ake + 2.839 Akzo = 0.13862 Akg3 


a ene rae cies ‘the aise: maximum error (p 
also the true maximum error, namely Az? = +1.129. ine 
Maximum Negative Ax’. As previously shown, the expression 
the primary error in 2° was 


; : 2Aku + 3Ak.2 — 0.05Ahi3 
zm ba = marty + 9(2Aka + 3Ake. — 0.05Ake3) 
+ 21(2Aks ae 3Ak39 a 0. 05Aks33) 


which, in order to yield its maximum negative value, requires the 
a omine signs: a 


ss ees Peet oad 
abe OS _ 

rr cats ee 
The primary errors obtained by the use of these signs are 6x (max.) 
= —1.050, dy = +0.169, and 62 = +0.093. 1 
The modified values of the unknowns will be approximately as — 


follows: 
xz’ = 2.000 — 1.050 = 0.950 


y’ = 3.000 + 0.169 = 3.169 
z’ = —0.050 + 0.093 = 0.043 


According to these indications, a change in sign occurs in z while x _ 
: is changing from z° to x’. Therefore, in solving for a tentative, or — 
apparent, maximum error, the signs given in the foregoing as being — 
required for obtaining the maximum primary error will be modified 


by reversing the signs of those Ak’s which are multiplied by this 3 
unknown (z), as follows: ’ 

- 

“hey / 

eo) oe oe : 

+ + + 

~ The revised equations are then 
2.20’ + 3.3y/ + 3.32’ = 12.85 ; 

l.la’ — 4.5y’ + 2.22’ = —13.10 

Liz! + 2.9y' — 0.97 =-_ 8.05 | 


| ee eG 


: Oreo or 1 (10: zat 1. 1( 22. 2.11): = 33.836 A, 


a 10.23)-++8.05( 22.11)]= 1.000, 
[12.85( 3.41) —13.10(—5.61)+8.05(— 1.21)]= 3.179 


gl: 85( 7.37)—13.10(—1. 21) +8.05(—13.53)] =0.0484 


The apparent maximum error in 2° is, Bacio: 1.000 — 2.000 
= —1.000 


ie 


1 
‘ Test: 62” = 33.836. (—0.79e; a 10.23¢, ea 22.11e;) 


where €; = 1.000AKia + 3.179Aki2 + 0.0484Aki3 


Since Ci, has changed sign, being now —0.79 as compared with 
+41.00 in the original equations, the signs of the Ak’s in ej should be 
reversed, giving the following combination of signs to be tried next: 


+ + + 
+ + + 


The new revised equations are: 


1.82" + 2.7y" + 2.72" 12.85 
Lie” — 4.5y" + 2.22" = —13.10 


ll 


l.1z’’ + 2.2y’’ — 0.92” = 8.05 
~ whence 
2D = 1.8 (—0.79)+ 1.1 ( 8.387)+1.1 ( 18.09) =27.684 
ge!’ =—_[12.85(—0.79) —13.10( 8.37)+8.05( 18.09)]= 0.933 
27.684 
1 
es 9g 12-85( 3.41) —13.10(—4.59) +8.05(—0.99)] = 3.467 


gi = : [12.85( 7.37) —13.10(—0.99) +8.05(—11.07)]= 0.670 
27.684 | | 


. 


kis 
: * 


‘ta 
Ak; in the expression for 62’’ 
Ky a 


It is now seen that every term 


ie 


a is negative for the signs used, namely, 


+ 4+ 4 


| + + 4+ = 
and therefore the true maximum error (negative) is Ar? = 0.933 
— 2.000 = —1.067. i, 


Example No. 10. Find the true maximum positive error in y?, 
using the data of Example No. 5. 
‘Solution: In Section 12 the expressions for the primary errors 
were found to be as follows: 


OF, 
Ey 
iz ee ut Ox |o 
‘ oe SS Pilar 
2 
r oe — a 
7 aF; 4 
| ae 4 
=. -and bx = ——_— ae \ 
ao —J OF, Ne 
E, 
OY |o 
OF, OF, 
where jae OR le 8Y Io 
OF, OF 4 
0% |o Oy Ne 
OF, 
= 2kex + ke = 12 
Ge |i 
OF, 


= 2(k3)2a° + kekey? = 28.5 


co} 
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Pt oF; 
/ oy ie Rs mae 
OF, 
ay [en hee? + Mie? = 31 


Ey => Ge yrAki + x°Ake + y°Aks ss bal 


and Ey = [2k3(x°)? + k3x°y?| Aks + [2k8(y°)2 + kay] Ak; = —20.45. 
The initial values of the unknowns were x° = 5, y° = 2; and the 
calculated primary errors were 6x = —0.144 and 6y (max.) = +0.792. 
The apparent maximum error will be calculated by using the 
same signs with the Ak’s as were required in obtaining the maximum 
positive primary error in y°, the modified values of the arbitrary 
constants being as follows: 
t= 1 1 = 11 
kg =2-+ 2 = 2.2 
ky = —34+.3 = —2.7 
ki = 1.5 — 15 = 1.35 


kg =2—-—.2=1.8 
In accordance with Equations (25) and (26), Section 12, the revised 


| equations are 
] 1.1(2’)? + 2.22’ — 2.7y' = 29 (40) 


1.8225(2')? + 2.432’y’ + 3.24(y’)? = 102.25 (41) 


In order to obtain the solution of these equations, the second 
method of Section 16 will be used. Approximate values of the 
_ unknowns, as found from the primary errors are 4.856 (for 2’) and 
2.792 (for y’). Using these values in the left-hand side of Equations 
(40) and (41), the following equations are obtained: 


11x? + 2.2% — 2.7y = 29.0836 
1.82252? + 2.43xry + 3.24y? = 101.1782 


By varying the right-hand terms of these equations to make them 
agree with Equations (40) and (41), and calculating the primary 
errors in the assumed values of the unknowns, the solution of the 
revised equations may be closely approximated. The required varia- 
tions in the right-hand terms are: Av; = —0.0836 and Av, = 1.0718. 


= Qkiz + ki, = 12.8832 


= 2(ki)2x + kikiy = 24.48468 


=k = -2.7 


oF. ; 
and : = = 2(ki)2y + kikix = 29.89224 
y 
The value of the Jacobian is therefore: 
.8832 —2.7 
= fee | = 451.2163 
24.48468  29.89224 
1 — 0.0836 —2.7 
h ‘cs = —___—_ = +0.0009 
ne” * 451.2163 hed eH . 
. 
1 12.8832 —0.0836 = 
and vy = | = +0.0351 
451.2163 |24.48468 1.0718 


The corrected values of x’ and y’ are thus found to be 4.8569 and 
2.8271, respectively, and when used in Equations (40) and (41) yield — 
the following results :* 


1.1(a’)? + 2.20’ — 2.7y’ = 29.00035 
1.8225(2’)? + 2.432'y’ + 3.24(y’)? = 102.25354 


The apparent maximum error in y° is therefore 2.8271 — 2.00 = 
+0.8271. ‘ 

In order to determine whether or not this is the true maximum 
error it is necessary to investigate the signs of the partial derivatives 


f 


*If greater precision in the solution of the revised equations is desired, a final set of corrections 
may be obtained by again varying the right-hand terms, putting An = —0.00035 and Avy = —0.00354. 


Pe SS ee 


—— 


weer et aes scehelig y ; Ne ih 1 
IN SIMULTANEOUS EQUATIONS 43 


RrOe abits, 
es ee using the modified values of the arbitrary constants 
5 


1 of the unknowns. The modified value of the Jacobian is as 
follows: 


2.22" + 2.2 ar 
3.6452’ + 2.43y’ 6.48y’ + 2.432’ 


U = 


- = 3 
124.57 30.12 


It is now seen that no changes have occurred in the signs of the 
minors of the Jacobian. Furthermore, since none of the unknowns 
have changed sign, it is evident that all of the error-coefficients have 
retained their original sign. A further increase in Ay can only be 
obtained by increasing the Ak’s beyond their assumed maximum 
values, and hence the true maximum error in y° is +0.8271. 


V. Errors IN FUNCTIONS OF THE UNKNOWNS 


18. Errors in Functions of the Unknowns.—The methods which 
have been used in the preceding sections, in connection with the 
calculation of errors in the unknowns, are applicable, with slight 
modifications, to the determination of errors in functions containing 


one or more of the unknowns 2... . Zn. 
Let 6 = h(a, 2%, ....; ki, ke, ....) be such a function; anid let 
=a, iy, 2-.. 3 1, Khe, --..).. The primary error in @° is. then 
_ given by 
00 06 00 00 
s0-( 621 = 622 + b, ear — Ak, Rr Pin +, 2) (42) 
4 O21 O22 Ok; , ° 
: In order to determine the maximum primary error in 6°, the initial 
_ values of the unknowns z?, 22, .... (found by solving the original 
| equations), and the expressions for 6x, 6%, .... are first obtained. 


_ These results are then substituted in Equation (42), and the terms 

are re-arranged so as to show each Ak with its proper coefficient. 

In this way, the expression for 69 is reduced to a linear function of 
the increments Ak; .... Ak», from which one can determine by in- 

_ spection the necessary signs which the Ak’s must have in order that 
the primary error in the given function shall be a maximum. 

The apparent maximum error may be computed by solving for 
a....x/,in the revised equations by one of the methods of Section 16, 
and then finding the corresponding value of 6’, whence A@? = 6’ — #°. 

If the true maximum error in @° is required, a criterion expression 
is found from the expression for 50 by replacing the initial values of 


eS 
Uy 
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‘a 
the arbitrary constants and of the unknowns by the corresponding 
modified values. The signs of the error-coefficients in the criterion — 


00 
expression, i. e., of ——.... 
: SE bells dk» 


described in Section 17 is followed in determining the signs and — 
magnitudes which the Ak’s must have in order to produce the true 
maximum error. 

The probable error in 6°, following the method described in Sec- 
tion 15, is obtainable from the terms in the expression for 60, after 
the latter has been reduced to the proper form for determining the 
maximum primary error. The probable error, \6°, is then found by 
taking the square root of the sum of the squares of the different 
terms in this expression. 

In order to illustrate the calculation of the maximum and probable 
errors in functions containing the unknowns of a system of simul- 
taneous equations, the following example, dealing with the bending 
moments in a continuous frame, will be solved. 


, are then examined, and the method 


Example No. 11. Find the maximum error in M,z, Fig. 1, due 
to possible variations of 10 per cent in the EJ/L values of the dif- 
ferent members; also find the probable error in M,, due to probable 
errors of +5 per cent in the EI/L values. 

The maximum moments will be computed for a uniform live load, 
the effects of dead load and side-sway being neglected. 


The relative span lengths and EI/L values are shown in the figure. 
Let M, denote the fixed-end moment in the middle span, and let 
k = EI/L. Also denote the rotations of joints X, Y, Z and W by 


' ee Sees 

Seo ae ee 
ae fag) tet = 
7  Keyt + 2(Key + kyy + kys)y + kysz 


oe 


i 


for which the following solution is obtained: 


asi 


ry ye ; 
= -M, (45) 
Kewd it 2(Kew + Keer) = 0 


(46) 


\ Vhen the Rie Priel k-values shown i in Fig. 1 are substituted, 


e equations become 
: bao + 2ye = 4M | (47) 
; 2x° + 14y° + 42° = —3M2 (48) 
4y° + 162° + 3? = —M; nc) 
32°+7we = 0 (50) 


5 2 0 0 
a “44 4 0 
Dim = 6238 
0 A 6 3 
0 0 3 7 
ees iy 
Mn if, 13814: 4) 0 
r= 
ie Ase 1(56) +0(—24)] =0.9429 M. 
Gi die Cay ows ie ‘ 
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Bae 23 Ol ectege [4(— 206) —3(515) 
"yea ead (eae ~1(—140)+0(60)] = —0.3573 Mp 
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*See Bulletin No. 108, Engineering Experiment Station, Unt of Ill., ‘‘Analysis of Statically 


Beane Structures by the Slope Deflection Method,” by 
amillo Weiss. d 


. M. Wilson, F. E. Richart, and 
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The function being studied is M. yay AS defined by the equatic n 
My = ky(2y + 2) — Ma 61) ) 


The initial value of this quantity is as follows: 


M3, (max.) = 4(—0.6854M2) — M. = —3.742M2 


= ang The expressions for the primary errors in the unknowns, as give n 


; by Equation (11b), are as follows: 
1 = 
» ~ 
1 ee 
02 may, 56e, —/ 140e, + 462e3 oe 198¢e.) (54) 


1 4 
dw = ae (— 24e,+ 60e. — 198e; + 976¢e,) . a 


From Equations (43) to (46), inclusive, the following values for 
é:.... @, are found: 


é: = 22°Akzy + (24° + y°) Akay 


= (2° + 2y°) Alay + 2y°Aleyy + (2y° + 2°)Akye 


rg 
th tithe) ras eee BE Te 


= (y° + 22°) Aky, + 22°Akeer + (22° + w°) Aken 
= (2° + 2w°)Akaw + 2w?Akwws 


wre 


ee as 


q Ra iies fee: 
= (1.886Akze + 1.528Akz,) Mo : (56) 

2 = (—0.715Akyy + 0.228Ak2z, — 0.685Aky:)M, (57) 

= (0.0584Ak.. — 0.299Ak,. + 0.046Akzn) Me (58) 

es = (—0.025Ak ww + 0.0042Ak.,) Mo (59) 


From Equations (42) and (51), the primary error in M®, is 
SMyz = ky.(2dy + 6z) + (2y? + 2°)Aky, 
or, since ky, = 4, and 2y° + 2° = —0.6854M, 
5M, = 8dy + 462 — 0.6854 Me Aky. (60) 


The value of the first two terms in Equation (60), as given by 
Equations (53) and (54), is as follows: 


8dy + 462 = 0.2285 e,; — 0.571 e, — 0.1167 e3 + 0.0501 e4 


and by using the values given in Equations (56) to (59), inclusive, 
_ and substituting in Equation (60), one obtains the following reduced 
expression for the primary error in Mj, in terms of the various Ak’s: 
6M, — M, (0.4381Ak.2 +0.408Akyy — 0.0068Ak,,2 me! 0.0012AK yy 
+0.220Ak,,—0.258Ak,,—0.0052Ak.) (61) 


From this expression it is seen that the specified variations in the 
k-values must have the following signs in order to produce the 
maximum primary error in M;,: 


Akze = —0.05, Akwu = +0.05, Akys = +0.4 
Akyy = —0.1, Aki = —0.2, Ake = +0.3 
Ake = +0.1, 


Upon substituting the foregoing values into Equation (61), the maxi- 
mum primary error in M2, is found to be —0.212Mz, which is equiva- 


lent to an error of 5.7 per cent. 
In order to determine the apparent maximum error, the revised 


equations are written and solved, namely, 
4.52 +1.8y’ = 4M, 
1.82’ + 14.2y’ + 4.42’ = —3M, 
4.4y’ + 17.62’ + 3.3w’ = —Mz 
3.32’ + 7.7w’ = 0 


= ——_- 


be ag 3 “hl 
n ' aa ia 
ee, _ 3(—224.334) — 1(60.984)] = 1.0298M@, 


vad 4, ‘em a) so 1. ae 224. 334) = = 6889, 2318 


eo (y! = ane — 3(560.835) — 1(—152.46)] = —0.3523M2 
a my ay 7 D’ aa 


= =" {4(60.984) — 3(—152.46) — 1(467.082)] = 0.0340M2 


w' = SHs(—26.186) — 3(65.35) — 1(—200.178)] = —0.0146M2 
. 
and M}, = ky.(2y’ + 2’) — Me 


= 4.4(—0.6707M2) — Mz = —3.951M; i 


The apparent maximum error is therefore 
AM3, = —3.951M2 — (—3.742M:2) = —0.209Mz 
Test: 5M), = ki,(2éy’ + 62’) + (2y’ + z’)Aky. 
a 8.85y’ + 4.462’ $37 0.6707 M2Ak,, 


; 
in eere (—224.384e; + 560.835¢, — 152.460; + 65.34e;) 


1 « 3 
5z' = an (60.984¢, — 152.46e, + 467.0826; — 200.178¢,) ~ 


where ef, = (2.060Akz2 + 1.708Aky) Me 
€, = (—0.702Akyy + 0.328Akzy — 0.668Aky.) Me 
= (0.068Ak.2 — 0.284Aky. + 0.053Ak.,)M2 
ek = (—0.029Akine + 0.005Aken) Me 


Also, 8.85y’ + 4.462’ = 0.2475e, — 0.619e4 — 0.1034e5 + 0.0476¢e; 
whence the reduced expression for 6M}, is as follows: 


5My.= (0.510 Ak +0.444Akyy —0.0070A Ke — 0.0014 A kyr 


,, : 
if ¥ . - 4 = 
ie ae Pees ae eee as ee ee ee 


+0.220Akzy—0.229AK yz — 0.0053 AK ew) Ma 


Since the separate terms in the foregoing expression will all be | 
negative when the Ak’s are given the signs which were used in com- 


apparent maximum error, the latter quantity is identical " 
1e true maximum error. 


Probable Error in My.. The probable error in M3,, found by sub- 
8 tuting the specified probable errors of 5 per cent toe the Ak’s in | 
the expression for the primary error, Equation (61), and taking the ~ 
hehe root of the sum of the squares of the resulting terms, is as 
ollows: 


[0.431(0.025)]? + [0.408(0.05)? + [0.0068(0.05) 
+ [0.0012(0.025)}? + [0.220(0.1)? | > 
+ [(0.258(0.2)]? + [0.0052(0.15)? ‘ 
= +0.065M, a 


VI. ConcLusions 


19. Comparison of Results—For purposes of comparison, the 
results of the preceding examples are tabulated below. 


Example Quantity Maximum Sadist M True Probable 
No. I ae | Primary aximum aximum 
" ae Error* Error* Error* aot 
4,7&9 x° = 2.00 +1.050 +1.129 +1.129 +0.298 
—1.050 —1.000 — 1.067 
5 &10 y? = 2.00 +0.792 +0.827 +0.827 +0.221 
1 it M3, = —3.742M,| —0.212M2| —0.209M.2| —0.209M.| +0.065Me 


*Assuming maximum errors of 10 per cent in the arbitrary constants. 
tAssuming probable errors of +5 per cent in the arbitrary constants. 


ao Vi pam 


20. Conclusions.—From the results of the foregoing investigation 
_ the following conclusions are offered: 

(1) The calculation of the maximum primary error and the 

probable error in a quantity governed by simultaneous equations is 
_ sufficiently simple to lend practical value to their determination. 

(2) When the errors in the arbitrary constants are in the neighbor- 
hood of 10 per cent, the maximum primary error will probably not 
differ from the true maximum error by appreciably more than 10 
per cent. 

(3) The apparent maximum error, although requiring more time 
for its determination than the maximum primary error, is a more 

satisfactory approximation to the true maximum error over a wide 


ye 


t of a Ak whose coefficient undergoes: 
changing from zm to 2m (or while a given function @ is ch: 
m 6° to 6’) may be considerable; prank one Spree ma: 


- certainty in estimating the Phagniiides of the errors in ‘ie ar bis 
_ trary constants. 

These conclusions are based upon a study of relatively strc ong 

systems of equations, and therefore should not be applied unreserved- 


ly to systems Spee ne the condition of linear dependence. * 


f 
.=. 


Le 
’ hone Elia A 


ee ae differentials of the functions Yeneaied ie 
ons (5) to (5n), Section 7, are as follows: 


= - i Gin Lm +. eran SHE mi ty 
+ 2dz,day + ary Aa + 2d2m0Ai1m + ery tac) - 2dx,dain = ()) (62) 


Shes 


= Ami0?xX1 +... . + Amm@tm +... + Amn@ Xn 
+2d2:dami+. ...+2d2mddmm+. ..-+2dtnddmn =0 (62m) 

Ofn = Anda +... . + Onmd?tm +... . + Anny 
+ 2daidan +... . + 2damddnm +.... + 2dtnddan = 0 (62n) 


24. 


= &2;, dx; 


oO 


= 62;, and dai; = Aa;;, 


By means of the relations 


oO 


the following equations are then obtained: 


aoa, +....+02,0%, +....+0%,8%, + >> dxjAa;=0 (63) 
i=l 


> 


aoa, +....+42,,0%r, +....+a%,02, + >) dajAa,; = 0 (63m) 
5 j=1 


| 02,8, +.... +.02,,0%, +....+02,8%, + >) bxjAan; = 0 (68n) 
j=1 


It will be noted that the determinant of the coefficients of the 
&zx’s in Equations (63) to (63n), inclusive, is identical with that of 
Equations (6) to (6n), in Section 7, and the terms preceded by the 
summation signs differ from the corresponding terms in Equations 


(6) to (6n), inclusive, only in having the unknowns 2... . Zn, in the 
latter, replaced by the primary errors, 6%... . 6%, ere ively, 


4 
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By letting de; = >> 62,;Aa.; (see Equation 7), the secondary error 


EM 
in 2?,, as found from Equations (63) to (63n), is then as follows: 


iy. <0 whe, Obst us 20S. 


(64a) | 


orgs Gr e.5 Oba ee ae 


Ops esie Dlg el eee 


or, using the co-factor notation, 


Whi n 
672m = ahiy DS C imde; (64b). 


an t=1 


In a similar way it can readily be shown that the errors 5°t,.... 
drm are given by 


ae ~. - Cimd?es (65) 
gta wa Se eanaie (66) 
—D i 
where ee, = so 6x ;Aa;; (67) 
and oN 416. = Y d4—ly Aa; (68) : 
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